Multiplexing biochemical signals 
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In this paper we show that living cells can multiplex biochemical signals, i.e. transmit multiple 
signals through the same signaling pathway simultaneously, and yet respond to them very specifically. 
We demonstrate how two binary input signals can be encoded in the concentration of a common 
signaling protein, which is then decoded such that each of the two output signals provides reliable 
information about one corresponding input. Under biologically relevant conditions the network can 
reach the maximum amount of information that can be transmitted, which is 2 bits. 



Cells continually have to respond to a myriad of sig- 
nals. One strategy for transmitting distinct stimuli is to 
use distinct signal transduction networks . It is, however, 
increasingly recognized that components are often shared 
between pathways [l|. Moreover, cells can transmit dif- 
ferent signals through one and the same pathway, and 
yet respond to them specifically. In rat cells, for instance, 
neuronal growth factor and epidermal growth factor stim- 
uli are transmitted through the same MAPK pathway, 
yet give rise to different cell fates, differentiation and 
proliferation respectively Q . These observations suggest 
that cells are able to transmit multiple messages through 
the same signal transduction network, just as many tele- 
phone calls can be transmitted via a single wire. Indeed, 
the intriguing question that arises is whether biochem- 
ical networks, like electronic circuits, can multiplex sig- 
nals: can multiple input signals be combined (encoded) 
simultaneously in the dynamics of a common signalling 
pathway, which are then decoded such that cells can re- 
spond specifically to each signal (see Fig. Q}? 

The question of how cells can transduce multiple sig- 
nals via pathways that share components is a key ques- 
tion in biology, since sharing components may lead to 
unwanted crosstalk between the different signals: from 
the perspective of one signal, the presence of additional 
signals constitutes noise. In recent years, several mech- 
anisms for ensuring signaling specificity have been pro- 
posed. One is spatial insulation, where the shared com- 
ponents are incorporated into distinct macromolecular 
complexes on scaffold proteins [l|. Other proposals are 
based on the temporal dynamics of the system, such as 
cross-pathway inhibition and kinetic insulation @. 
However, these studies only considered scenarios in which 
the system is stimulated with one signal at the time. 
Rensing and Ruoff studied what happens when two or 
three MAPK pathways that share components are stim- 
ulated simultaneously jjj], but found that one pathway 
tends to dominate the response, suggesting that multiple 
messages cannot be transmitted simultaneously. Here we 
demonstrate that cells can truly multiplex signals: we 
show that they can transmit at least two signals simul- 
taneously through a common pathway, and yet respond 



specifically to each of them. 

We first have to understand how multiple signals can 
be encoded in the dynamics of a signaling pathway. Cells 
employ a number of coding strategies for transducing sig- 
nals. One is to encode stimuli in the temporal dynamics, 
such as the duration Q or frequency [§], of an intracel- 
lular signal. In principle, these coding strategies could 
be used to multiplex signals. Here, we consider what is 
arguably the simplest and most generic coding strategy 
cells could choose, namely one in which the signals are 
encoded in the concentrations of the signaling proteins. 
We will call this strategy AM multiplexing. 

We will consider the biochemical network shown in 
Fig. UK- It consists of N input species Si, . . . , Sat with 
copy numbers Si, . . . , Sn, a signal transduction path- 
way V consisting of M species Vi, . . . , Vm, and N out- 
put species Xi, . . . , Xjy. The copy number of each input 
species Sj can be in one of K states, = 0, . . . , K — 1, 
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FIG. 1: (a) Biochemical multiplexing: N different signals are 
encoded in the state of a common pathway V, which is then 
decoded such that each output species Xi provides reliable 
information about the corresponding input Si. (b) Multiplex- 
ing is a mapping problem. The states of two inputs Si and 
S2 are mapped onto the concentration of V, which is then 
mapped onto states of the output species Xi and X2; we re- 
quire that the two lowest (highest) levels of Xi correspond to 
the lowest (highest) level of Si; the dashed arrow denotes a 
mapping that violates this requirement; levels of V and Xi 
are colored according to input pattern s = (si, sa)- (c) The 3 
unique mappings of s to v; in panel (b) mapping C is shown. 
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which are labelled in order of increasing copy number, 
S<l 0) < Sf } < ■•■ < Sf^ 1 ). The input pattern is de- 
noted by the vector s = (si, . . . , Sjv)- Similarly, the copy 
number of each output species Xj can be in one of L 
states x-i = 0, . . . , L — 1 ordered by increasing copy num- 
ber Xj, and the output pattern is denoted by the vector 
x = [xq, . . . , xjv)- A necessary condition for multiplex- 
ing is that the state space of V is large enough that it 
is possible to encode the total number of input patterns, 
K N , in V. 

We imagine that the N input signals are independent, 
and that the signal transduction network V replaces N 
independent signaling pathways. We therefore require 
that Xi should provide reliable information about the 
state Si, but not necessarily about Sj^; the N different 
input signals s simply have to be transduced to x, not 
necessarily integrated. In general, however, the state Xj 
will be a function of the states of all the input species: 
Xi = f(s). This reflects the fact that inevitably there is 
cross-talk between the different signals because they are 
transmitted via the same pathway. However, this cross- 
talk is not detrimental as long as it does not compromise 
the cell's ability to infer from Xi what Sj was. 

Another key point is that while the precise mapping 
from s to x may not be critical for the amount of informa- 
tion transmitted per se, this is likely to be important for 
whether or not this information can be exploited. Let's 
imagine that the system contains three input species, say 
three sugars, and that each of these can be in one of only 
two states, Sj = or 1, corresponding to the absence 
or presence of the sugar; let's further assume that Xj is 
an enzyme needed to consume sugar Sj. With 8 input 
patterns Xi can, in the absence of noise, take 8 values, 
identified as states Xj = 0, ... ,7. Now, it seems natural 
to demand that when the sugar Sj is absent (sj = 0), 
the copy number of enzyme Xj is low, while when Sj is 
present, the copy number of Xj is high; this means that 
the four lowest levels of Xi (xj = 0,1,2,3) should cor- 
respond to Si = 0, while the four highest levels of Xi 
should correspond to s, = 1. We therefore require that 
the mapping from s to x is such that the output states 
{xi} corresponding to input Sj = j are grouped into sets 
that are contiguous and either increase or decrease mono- 
tonically with j, for each signal i. This leads to a mono- 
tonic input-output relation between Si and Xi for each i. 
We call this requirement the multiplexing requirement. 

In the rest of the manuscript, we make these ideas con- 
crete for a network in steady state with two input species, 
Si and S2, each of which has either a low (s, = 0) or a 
high concentration (sj = 1). We take a signaling path- 
way V consisting of only one species, V. Multiplexing 
requires that, in the absence of noise, the four input pat- 
terns s can be mapped onto four distinct states of V, 
v = 0, . . . , 3, again labelled in order of increasing copy 
number. These four levels of V lead to four states for 
each of the two output species Xi and X 2 (Fig. QJ3). As 



explained above, we require that we can group these four 
states into two sets, called LOW and HIGH, such that 
the LOW set, containing Xj = 0, 1, corresponds to Sj = 
and the HIGH set, containing Xj = 2, 3, corresponds to 
Si = 1 (or vice versa, leading to an inverse input-output 
relation). To elucidate which mechanisms make it pos- 
sible to multiplex Si and S2, we note that there exists 
different ways of mapping s to v, but, as we will explain 
shortly, not all of these mappings can be decoded into x 
in a manner that satisfies the multiplexing requirement. 
We therefore first address the question which combina- 
tions of mapping from s to v and decoding from v to 
x fulfill the multiplexing requirement, and then we will 
discuss what encoding mechanisms actually allow for the 
required mapping from s to u. 

Due to the symmetry in the problem, there are 3 
unique ways of mapping the four input patterns s to v 
(Fig. HP). To determine whether there exists a scheme 
for decoding the signals from v to x that satisfies the 
multiplexing requirement, we examine for each mapping 
all possible network topologies between V, Xi and X 2 , 
except those that involve autoregulation or mutual re- 
pression/activation since these may lead to bistability. 
In particular, we allow not only for activation and re- 
pression of Xi and X 2 by V, but also for activation and 
repression of X 2 by Xi, leading to feedforward loops, a 
common motif in signal transduction pathways and gene 
networks Q. In the deterministic mean- field limit the 
steady-state values of X\ and X 2 are thus given by 

X 1 = kxf{V;K ai n a )/n, (1) 
A 2 = k 2 f{V;K ,np) X f(X 1 ;K 1 ,n y )/^, (2) 

where k is the maximum activation/production rate, 
fi is the degradation/deactivation rate, and each reg- 
ulation function is either an activating or repress- 
ing Hill function, f{V;K,n) = V n /(V n + K n ) or 
f{V;K,n) = K n /(V n + K n ). The multiplication in 
Eq. [2] indicates that we assume that at X 2 , X\ and 
V are integrated according to AND logic [9j. We per- 
formed extensive sampling of the space of parameters 
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for each of the mappings 



ki,k 2 ,K a ,n a ,Kp,np, 
in Fig. fflC. 

Only for mapping C do we find decoding schemes that 
satisfy the multiplexing requirement 16j . Interestingly, 



all valid decoding networks are incoherent feedforward 
loops [9(. Figure [2] illustrates the principle for one such 
motif. Panel B shows for each of the four input patterns 
s the copy number V together with the threshold copy 
numbers, K a , Kp and K 7 , while panels C and D show 
X\ and X 2 respectively as a function of V. Xi(V) is 
a simple activation curve with activation threshold K a . 
In contrast, X 2 (V) starts low and rises around Kp, but 
then decreases again due to repression by Xi. This non- 
monotonicity, which is a result of the incoherent charac- 
ter of the feedforward loop, is critical, since this makes 
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FIG. 2: Decoding V using an incoherent feedforward loop, 
(a) Network architecture, (b) The values of V corresponding 
to the four input patterns s for mapping C (see Fig. QTJ), with 
thresholds K a ,K ft ,K 1 (see Eqs. [T|2). (c) Xl(V); (d) X 2 {V). 
The non-monotonicity of X2(V) swaps the states correspond- 
ing to (1, 1) and (0, 1) in the mapping from v to 



it possible to swap the order of the states corresponding 
to s = (1,1) and (1,0) in the mapping from v to x 2 . The 
key parameters are the activation/repression thresholds 
K, since they determine where in the state space of V 
the outputs switch between high and low levels. The pre- 
cise values of k and n are of less importance, although n 
should not become so large that X\(V) becomes Boolean: 
it is critical that Xi, which needs to be activated by V 
around K a to transmit S\ 1 is not fully activated at K a : 
to multiplex 5*2 , X\ should reach the threshold K 1 for re- 
pressing X2 only when V has become significantly larger 
than K a . Indeed, if X\ can only take two states, then 
only three states of V could be decoded, and not the re- 
quired four. AM multiplexing thus relies on the fact that 
signals can be encoded over a range of concentrations. 

We can now also understand why mappings A and 
B are difficult to decode: they would require an input- 
output relation between X2 and V that rises more than 
once. This is difficult to achieve in a feedforward loop 
without mutual repression or activation. 

The above analysis shows that it is possible to decode 
multiple signals simultaneously, provided that the input 
s can be encoded in v according to mapping C. The next 
question is how these mappings, which correspond to 
particular input-output relations V(Si , S 2 ). can be gen- 
erated. Experiments [l(| and modelling [ill 12 1 have 
shown that transcriptional regulation can be very sophis- 
ticated, allowing for complex logical operations 12|. We 
indeed find that a simple scheme for transcriptional reg- 
ulation based on the mechanism of 'regulated recruit- 
ment' [ll| can generate the required input-output rela- 
tion V{S\,S 2 ), where Si and S2 are now transcription 
factors that regulate the expression of the protein V. In 
this scheme, Si and S2 independently activate gene ex- 
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FIG. 3: (a) An input-output relation V(Si,S2) consistent 
with mapping C (Fig. HP); symbols correspond to 

states in Fig. [2] (b) Architecture of the promoter, (c) A 
feedforward loop that can generate mapping C. 



pression by binding next to the core promoter, thus re- 
cruiting the RNA polymerase (RNAp), while Si and S2 
together repress gene expression by cooperative binding 
to the core promoter, thereby blocking the binding of 
RNAp (see Fig.©. This yields 



V(S 1 ,S 2 ) 



(/?///) g p (l + coqi + uq 2 ) 



l + q[+q 2 + uj'q[q 2 + q p (l + UJqi + coq 2 ) ' 

(3) 

where j3 is the maximum expression rate and /j, is the 
degradation rate of V, q p = c p /K p is the concentration of 
RNAp c p scaled with its dissociation constant K p , q\ = 
Si/Kt, q 2 = S 2 /K 2l q[ = Sx/Ki, and q' 2 = S 2 /K 2 , 
where Ki and K[ are the dissociation constants for the 
binding of Si to the promoter sites where the RNAp is 
recruited or blocked, respectively; uj and u>' are factors 
reflecting coop erative interactions between the respective 
molecules [llLll2l] . We thus conclude that gene regulation 
networks have the capacity to multiplex signals. 

While it is clear that signaling pathways often share 
common components [ILQ, the logic of signal integration 
in these pathways has been characterized in much less de- 
tail than for gene regulatory networks. It is conceivable 
that the desired input-output function V(Si,S 2 ) could 
be implemented at the level of a single protein V, us- 
ing competitive and/or cooperative binding between the 
three molecules Si, S2, V. Alternatively, the required 
encoding could also be implemented at a higher level of 
network interactions. For instance, a network in which Si 
and S2 regulate V via two additional components, Qi and 
Q 2: in an incoherent feedforward loop (Fig. [3]C), could 
achieve the required encoding V(S±, S 2 ). In essence, the 
feedforward loop between Qi, Q 2 and V can be used to 
control the ordering of V in the encoding process, just as 
the feedforward loop between V, Xi and X 2 can be used 
to regulate the ordering of X 2 in the decoding step. Since 
feedforward loops are common motifs in signal transduc- 
tion pathways [9( , we argue that multiplexing can also be 
implemented in these networks. 

The analysis above shows that in principle biochemical 
networks can multiplex signals in the mean- field, deter- 
ministic limit. However, there remains the question of 
whether signals can be multiplexed reliably in the pres- 
ence of inevitable biochemical noise. To address this, we 
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FIG. 4: The transmitted information I as a function of y" max 
for four different values of X max . / can reach the maximum 
value of 2 bits provided y max and X max are large enough. 

estimate a lower bound on the information about two 
binary signals Si and S2 that are transmitted through 
the network studied above (Eqs. [T]|3]). We define the to- 
tal information / = I(Si,Xi) + I(S2, X2) as the sum of 
the mutual information for each of the individual signals, 
US l ,X l ) = ^2 x .J2 Si P( x h s i) io s\p(xi,Si)/p(xi)p(si)] 
jl3l |. where p(si) and p(xi) are respectively the probabil- 
ities of Si being in state s$ and Xi being in state Xi , and 
p(xi, Si) is the joint probability of input Sj and output x^. 
Note that in the presence of noise Xi is not limited to 4 
states but can in principle take any value. This definition 
of / makes it straightforward to directly compare the per- 
formance of this network with that of two independent 
pathways. If each of the two input states for each Si is 
equally likely then the maximum value of I(Si,Xi) is 1 
bit for each transmitted signal i; the maximum value of 
I is thus 2 bits. 

To maximize the lower bound on / we optimize the net- 
work parameters using a simulated-annealing algorithm; 
we have verified that the final results are robust by vary- 
ing the initial conditions, and by also using an evolution- 
ary algorithm. We fix the degradation rate of all proteins 
to be /i = lhr -1 and vary n a , np and n 7 between 1 and 
4. Values of k\, are set such that the maximum mean 
value of each Xi is AT max ; similarly, (3 is set such that the 
maximum mean value of V is U max . x max and y max are 
varied systematically (see Fig. 2]). The threshold param- 
eters K a , Kp and K 1 are varied over the range [0, y max ] 
or [0,X max ] as appropriate. We vary q p , qi from 10~ 2 
to 10 2 and ui, uj 1 between 1 and 10. For each parameter 
set we compute p{xi,Si) using the linear-noise approx- 
imation [14j . Its accuracy was verified by performing 
Gillespie simulations of the optimized networks [151 ] . 

Figure |4] shows that below a threshold copy number 
ym&x QQ j-jjg totaj information is low regardless of 
X max because four distinct states of V cannot be gen- 
erated. Above V^, max , for large X max the information I 
reaches 2 bits, the maximum information about the two 
signals S% and S2 that could be transmitted via two inde- 



pendent channels. For lower values of X max , / saturates 
at a value lower than 2 bits, limited by the intrinsic noise 
in the production and decay of Xj. Importantly, I reaches 
2 bits for V" max w X max ps 500, which is well within the 
range of typical protein copy numbers inside living cells. 
This shows that biochemical networks can multiplex two 
signals reliably in the presence of biochemical noise under 
biologically relevant conditions. 

In summary, our results suggest that cells can trans- 
mit at least two binary signals through one and the same 
pathway, and yet respond specifically and reliably to each 
of them. The proposed mechanism for biochemical mul- 
tiplexing is based on swapping the order of states during 
the encoding and decoding steps using incoherent feed- 
forward loops. It is clear that the principle is generic, 
and could be implemented in signal transduction path- 
ways and gene networks - indeed incoherent feedforward 
loops are commonly found in these networks [9]. Our 
predictions could be tested experimentally by simultane- 
ously stimulating two MAPK pathways that share com- 
ponents although perhaps a more controlled exper- 
iment would be one using synthetic gene networks. In 
future work, we will address how more than two input 
signals can be transduced simultaneously, and how cells 
can multiplex signals by encoding them into the temporal 
dynamics of the signaling pathway. 
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